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Nitrate is an inliibitor of sulfate-reducing bacteria (SRB). In petroleum production sites, 
amendments of nitrate and nitrite are used to prevent SRB production of sulfide that 
causes souring of oil wells. A better understanding of nitrate stress responses in 
the model SRB, Desulfovibrio vulgaris Hildenborough and Desulfovibrio alaskensis G20, 
will strengthen predictions of environmental outcomes of nitrate application. Nitrate 
inhibition of SRB has historically been considered to result from the generation of small 
amounts of nitrite, to which SRB are quite sensitive. Here we explored the possibility 
that nitrate might inhibit SRB by a mechanism other than through nitrite inhibition. 
We found that nitrate-stressed D. vulgaris cultures grown in lactate-sulfate conditions 
eventually grew in the presence of high concentrations of nitrate, and their resistance 
continued through several subcultures. Nitrate consumption was not detected over 
the course of the experiment, suggesting adaptation to nitrate. With high-throughput 
genetic approaches employing TnLE-seq for D. vulgaris and a pooled mutant library of 
D. alaskensis, we determined the fitness of many transposon mutants of both organisms 
in nitrate stress conditions. We found that several mutants, including homologs present 
in both strains, had a greatly increased ability to grow in the presence of nitrate but not 
nitrite. The mutated genes conferring nitrate resistance included the gene encoding the 
putative Rex transcriptional regulator (DVU0916/Dde_2702), as well as a cluster of genes 
(DVU0251-DVU0245/Dde_0597-Dde_0605) that is poorly annotated. Follow-up studies 
with individual D. vulgaris transposon and deletion mutants confirmed high-throughput 
results. We conclude that, in D. vulgaris and D. alaskensis, nitrate resistance in wild-type 
cultures is likely conferred by spontaneous mutations. Furthermore, the mechanisms that 
confer nitrate resistance may be different from those that confer nitrite resistance. 
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INTRODUCTION 

Sulfate-reducing bacteria (SRB) are environmentally and indus- 
trially significant microorganisms that use sulfate as a terminal 
electron acceptor in anaerobic respiration. These anaerobes pro- 
duce sulfide as the end product of sulfate respiration (Postgate, 
1984). Sulfide is toxic to most organisms (Caffrey and Voordouw, 
2010), and its production causes oil souring in the petroleum 
industry (Ligthelm et al., 1991; Sunde et al, 1993). Despite the 
undesirable features of this metabolic end product, SRB have been 
exploited in studies of heavy metal bioremediation (Jiang and 
Fan, 2008; Martins et al, 2009) because of the ability of sulfide to 
form insoluble complexes with heavy metals (Jalali and Baldwin, 
2000). SRB also precipitate heavy metals by directly changing the 
metal redox state to a less soluble form (Lovley et al., 1993a,b; 
Lloyd et al., 1999; Chardin et al, 2003). The metabolism of SRB 
is studied, therefore, to understand how to minimize the detri- 
mental economic effects of these bacteria and to maximize their 
positive metabolic traits. 



Such traits have been studied extensively in Desulfovibrio vul- 
garis Hildenborough (DvH), a model SRB with a sequenced 
genome (Heidelberg et al, 2004). DvH has been examined under 
a variety of stress conditions, including elevated nitrite (He et al., 
2006; Bender et al., 2007) or nitrate concentrations (Redding 
et al, 2006; He et al, 2010a), heat shock (Chhabra et al, 2006), 
high salt (Mukhopadhyay et al., 2006; He et al, 2010b), oxygen 
(Mukhopadhyay et al, 2007), or electron donor depletion (Clark 
et al., 2006). The data obtained help in the prediction of responses 
of SRB in heavy metal-contaminated sites, which also contain 
many chemicals that inhibit these bacteria. For example, nitrate 
concentrations can be greater than 100 mM at US nuclear weapon 
complexes overseen by the Department of Energy (Green et al., 
2012), and these waste sites are also contaminated with heavy 
metals (Riley and Zachara, 1992). High nitrate inhibits the growth 
of DvH (He et al., 2010a). Although some SRB can also use nitrate 
as a terminal electron acceptor (McCready et al., 1983), nitrate is 
successfully used by the petroleum industry to control the growth 
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of SRB and the oil souring that their sulfide production causes 
(Sunde and Torsvik, 2005). The mechanism of nitrate inhibition 
of SRB is still unclear. In the environment, at least part of the inhi- 
bition by nitrate is indirect: nitrate-reducing, sulfide-oxidizing 
bacteria produce nitrite that is toxic to SRB at much lower con- 
centrations than is nitrate (Haveman et al., 2005; He et al, 2006). 
Furthermore, in oil wells, heterotrophic nitrate-reducing bacteria 
can compete with SRB for volatile fatty acid electron donors, fur- 
ther reducing the production of sulfide (Grigoryan et al., 2008). 
However, nitrate is also inhibitory to DvH in the absence of 
nitrate-reducing bacteria (Redding et al., 2006; He et al., 2010a). 

It has been suggested that this pure culture nitrate inhibition 
is also a result of nitrite stress, since DvH itself may produce 
small amounts of nitrite from non-specific reduction of nitrate 
(Wall et al., 2007). In addition, high concentrations of nitrate 
could potentially induce a non-specific osmotic shock response 
in the bacteria (Wall et al., 2007). However, microarray data reveal 
few common gene expression changes among nitrate, nitrite, and 
sodium chloride stress conditions (He et al., 2010a). He et al. sug- 
gested that unique nitrate stress responses might account for these 
discrepancies (He et al., 2010a). Understanding the mechanism of 
nitrate inhibition of DvH and the genes involved in the nitrate 
stress response should facilitate the prediction and monitoring 
of the effectiveness of bioremediation strategies that employ SRB 
(Hazen and Stahl, 2006). 

Past studies of the mechanisms of nitrate stress responses in 
DvH have relied primarily on transcript analyses (He et al., 2010a) 
and protein determination (Redding et al., 2006) techniques. 
However, mutant analysis is a more reliable method of determin- 
ing gene essentiality in a particular stress condition (Price et al., 
2013). Fitness profiling of many mutants en masse is a high- 
throughput approach complementary to classical genetic tech- 
niques that has allowed rapid annotation of genes (Deutschbauer 
et al., 2011). In this study, we used random transposon mutant 
fitness profiling in two completely sequenced (Heidelberg et al., 
2004; Hauser et al., 2011) model SRB, Desulfovibrio alaskensis 
G20 ("G20," formerly called Desulfovibrio desulfuricans G20) and 
DvH, to probe the molecular mechanisms of their nitrate stress 
responses. Because 58% of G20 genes (1954/3371) are also present 
in DvH, which has 3503 genes (http://www.microbesonline.org/), 
we predicted that the strains would have similar nitrate stress 
responses. Therefore, pools of DvH and G20 transposon mutants 
with mutations saturating the non-essential genes under permis- 
sive growth conditions were subjected to high concentrations of 
nitrate. Illumina sequencing or microarrays were used to locate 
the transposons in mutants surviving the nitrate exposure and, 
by comparison with mutants not exposed to stress, to iden- 
tify genes essential for survival in high nitrate. Generally, those 
mutants lost from the stress treatment represent genes whose 
functions are needed for stress survival. From the fitness profiling 
reported here, surprisingly we identified mutants with dramati- 
cally increased fitness in nitrate stress conditions that we further 
analyzed in pure cultures. However, the same mutations did 
not confer resistance to nitrite. These results confirmed the pre- 
dicted existence of unique nitrate-resistance mechanisms (He 
et al., 2010a) and suggested that environmental models of nitrate 
inhibition require expansion. 



MATERIALS AND METHODS 
STRAINS AND MEDIA 

The strains used in this study are listed in Table 1. All DvH and 
G20 strains were grown in defined MOLS4 medium [MO Basal 
Salts (Zane et al, 2010) with 60 mM sodium lactate and 30 mM 
sodium sulfate]. The medium used to grow DvH cultures was 
reduced with 1.2 mM sodium thioglycolate; whereas, the medium 
for G20 was reduced with 0.38 mM titanium citrate. DvH and 
G20 manipulations, including setup of growth kinetic studies, 
were done at about 25°C in an anaerobic growth chamber (Coy 
Laboratory Products, Inc., Grass Lake, MI) with an atmosphere 
of approximately 95% N2 and 5% H2. Optical densities (600 nm) 
were determined with a Genesys 20 spectrophotometer (Thermo 
Scientific, Waltham, MA). 

GROWTH KINETICS 

Growth kinetic studies with WT G20 were carried out essentially 
as described below for the G20 fitness profiling, with the fol- 
lowing exceptions: No kanamycin was used with the WT cells, 
and each wild-type G20 freezer stock was pelleted to reduce car- 
ryover of glycerol used as the cryoprotectant before inoculation 
of starter cultures. For all DvH growth kinetics, 5 mL MOLS4 
cultures (with 1.2 mM sodium thioglycolate) were started by 
inoculation with the pelleted cells from a freezer stock. Geneticin 
(G418) sulfate (400|jLg/mL) or spectinomycin dihydrochloride 
pentahydrate (100(JLg/mL) were added to DvH cultures where 
indicated. Each condition tested was prepared as 14.5 mL of inoc- 
ulated culture plus 1 mL deionized water (for "no additions" 
controls) or inhibitory salts (sodium nitrate, sodium nitrite) dis- 
solved in deionized water. Aliquots of 5.1 mL from this 15.5 mL 
culture were grown as triplicates in 27-mL anaerobic tubes, each 
capped with a butyl rubber stopper and grown at 34°C. All G20 
and DvH inocula were grown to late exponential or stationary 
phase. 100 mM sodium nitrate was used for experiments with 
DvH, compared with 150 mM for G20 experiments, due to the 
greater sensitivity of DvH to nitrate. With a few exceptions, all 
growth kinetics experiments were repeated at least twice with 
triplicates in each experiment. Triplicate growth experiments for 
the DVU0250 transposon mutant and intergenic control transpo- 
son mutant experiments were done once. Triplicates of the A rex 
mutant in the presence of nitrite were also grown once. 

G20 TRANSPOSON MUTANT FITNESS STUDIES 

Fitness data were collected with two pools of G20 transposon 
insertion mutants (4069 and 4056 mutants, respectively) that will 
be described in more detail elsewhere (Kuehl et al., in revision). 
Briefly, 1 174 strains are present in both pools, leaving 6951 strains 
that are present only once in the complete library. 498 genes are 
represented only once in the library. 571 genes are represented 
twice in the library; that is, either a single strain is present in 
each of the pools or two different strains with a transposon inser- 
tion in the same gene are present in the library. 1341 genes are 
represented in the library three or more times. A total of 2410 
unique genes and 212 unique intergenic regions are represented. 
Thus, about 71% of G20 genes are represented in the library, 
providing excellent coverage of non-essential genes. Transposon 
insertions were mapped to the genome by a two-step arbitrary 
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Table 1 | Strains and plasmids used in this study. 



Strain or plasmid 



Genotype or relevant characteristics^ 



Sources and/or references 



a-Select (Silver Efficiency) 



D. alaskensis ^ 

G20 

D. vulgaris 

ATCC29579 
JW710 

JW3311 

GZ9685 

GZ12997 

GZ2640 

GZ12015 

GZ10694 

GZ2179 



PLASMIDS 

plVI0719 

pl\/10746 

pl\/lO9075 

pIVlOSSII 



pIVlOSSIS 
pRL27 



F" deoR endA^ recAl re/AI g/rA96 hsdR^7{r~, 
A(/acZYA-srgF)U169 *80/scZAM15 X." 



I supE44 thi-^ phoA 



Spontaneously nalidixic acid-resistant derivative of Desulfovibrio desulfuricans 
G100A lacking the endogenous cryptic 2.3-kb plasmid, pBGI 



pCR8/GW/T0P0 containing SRB replicon (pBGI); Sp^■ source of Sp^ and pUC ori 
fragment for marker exchange suicide plasmid construction 
upp in artificial operon with npt and Ap'"-pLIC ori; Pnpt-npt-upp; Km'; 5-FLI^; for 
marker exchange suicide plasmid construction 

pM0719 containing Pnpt for constitutive expression of complementation 
constructs: pBGI stable SRB replicon; Sp'" 

Sp'^ and pUC ori from pM0719 plus 1630 bp upstream and 1590 bp downstream 
DNA regions from DVU0916 {rex) flanking the artificial operon of Pnpt-npt-upp 
from pM0746; for marker exchange deletion mutagenesis; Sp'" and Km'^ 
pMO9075 with DVU0916 (rex) constitutively expressed from Pnpt 
Tn5-RL27 (Km'"-oriR6 K) delivery vector; for transposon mutagenesis of DvH 
strains 



Bioline 



Walletal., 1993 



Wild-type D. vulgaris Hildenborough (pDVI ); 5-FU= (Parent for GZ strains) ATCC 

WT Aupp (pDVI); 5-FU'" (used as "WT" control for DvH growth kinetics in this Keller et al., 2009) 

study; parent strain for deletions) 

JW710 ADVU0916::(npt upp); Km'; 5-FU^(Arex marker exchange) This study 

DVU0245-773::Tn5-RL27;insertion at bp 773/1110 for the gene; Km^ Wall laboratory 

DVU0246-111::Tn5-RL27;insertion at bp 111/2235 for the gene; Km^ Wall laboratory 

DVU0247-211::Tn5-RL27;insertion at bp 211/360 for the gene; Km^ Wall laboratory 

DVU0250-427::Tn5-RL27;insertion at bp 427/588 for the gene; Km^ Wall laboratory 

DVU0251-80::Tn5-RL27;insertion at bp 80/963 for the gene; Km^ Wall laboratory 

Genome position 658487::Tn5-RL27; insertion at intergenic region 327 bp Wall laboratory 
upstream of VIMSS209534, DVU0590; Km' (Control strain for transposon mutant 
growth kinetics) 



Keller etal., 2009 
Parks etal., 2013 
Keller etal., 2011, 2014 
This study 



This study 
Larsen et al., 2002 



^Km, kanamycin; Sp, spectinomycin: Ap, ampicillin: 5-FU, 5-fluorouracil; superscript "r" or "s," resistance or sensitivity. 



PGR as described previously (Oh et al, 2010). Each mutant has 
a "TagModule" that contains two different variable segments, an 
"uptag" and a "downtag"(Oh et al., 2010). Within each pool, each 
strain has a unique TagModule, so that the abundance of the 
TagModule is a proxy for the abundance of that strain. Only the 
uptags are amplified from the "uptag pool," Pool 1 (4069 strains) 
and only the downtags are amplified from the "downtag pool" 
Pool 2 (4056 strains). Amplified tags from both pools can be 
hybridized to the same array because only one tag (up or down) 
from a TagModule has been shown to be necessary for accurate 
quantification of strain abundance and there is no overlap of tags 
in the two pools (Oh et al, 2010; Deutschbauer et al, 2011). 
Each pool was grown overnight to late log phase (OD about 0.87) 
in about 10 mL MOLS4 medium amended with titanium citrate 
(0.38 mM) as reductant and kanamycin (800 (Jig/mL). 750 (jlL of 
each pool was added to 15 ml MOLS4 + 1650 |xL sterile MO 
Basal salts (Zane et al., 2010) or salt (sodium nitrate, sodium chlo- 
ride, etc.) dissolved in MO Basal Salts. Each amended medium 
plus mutants (17.4 mL) was aliquoted into three 27-mL anaerobic 



tubes, about 5.8 mL per tube, each capped with a butyl rub- 
ber stopper and grown at 34° C. When the cultures had reached 
stationary phase (OD > 1), 0.5 mL from each control or stress 
condition was pelleted and processed as described previously; that 
is, genomic DNA was extracted (Deutschbauer et al., 2011), and 
the uptags and downtags were PGR amplified, hybridized to an 
Affymetrbc 16K TAG4 array, and scanned (Pierce et al., 2007). 
The number of doublings of the population was estimated by 
using the doubling in ODsoo to approximate doubling of the cell 
population. 

Fitness data for G20 were analyzed as described for simi- 
lar experiments with Shewanella oneidensis MR-1 (Deutschbauer 
et al., 2011), with slight modifications (Price et al, 2013). Briefly, 
strain fitness = log2 (END/START), where those values ("END" 
and "START") are averages of the gene location-specific uptag 
and downtag log2 intensities. Mutant strains with low START val- 
ues were excluded, leaving measurements for 3726 strains in Pool 
1 and 3865 strains in Pool 2. Strain fitness was normalized across 
the genome so that the median was 0; this was done separately 
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for the two pools. Since a gene could be mutated at different sites, 
gene fitness was calculated as the average fitness of strains with 
mutations in a particular gene. Gene fitness was normalized to 
remove the effect of chromosomal position on gene fitness and 
to set the mode of fitness values to zero, as previously described 
(Price et al., 2013). One difference from the previously described 
protocol (Price et al., 2013) was that only one Affymetrix chip 
was used per experiment; up- and downtags were hybridized to 
the same array. The results of the "MOLS4 no stress" condition 
were similar to the LS4D controls described previously (Price 
et al, 2013), including similarly "sick" auxotrophs that would 
be expected in a defined medium when compared with lactate- 
sulfate medium containing yeast extract. Two quality metrics were 
used for each experiment. Strain correlation (Table SI) is the cor- 
relation of the strain fitness values for the same strains between 
the two pools. Operon correlation (Table SI) is the correlation of 
gene fitness values between adjacent genes predicted (http://www. 
microbesonline.org/) to be in the same operon. The low quality 
metrics for the sodium nitrate and potassium nitrate conditions 
(Table SI) reflect the predominance of a few strains in the cul- 
ture; essentially only the data for these few strains is biologically 
meaningful, while the majority of strains did not have the oppor- 
tunity to grow at all before the culture reached stationary phase. 
The complete data from these experiments are available at http:// 
www.microbesonline.org/. 

DvH TnLE-Seq FITNESS STUDIES 

This nitrate resistance study was one of five multiplexed TnLE-seq 
pools that were part of a protocol described elsewhere (Pels et al., 
2013). One advantage of this method is that individual mutants 
do not need to be isolated and confirmed, as in a catalogued 
library. They also do not need to be frozen en masse and recovered 
from the freezer. Rather, hundreds of thousands of unique trans- 
poson mutations are created by conjugation at the beginning of 
each experiment. The only experimental difference between this 
study and those published was the addition of 100 mM nitrate to 
the MOLS4 defined medium for growth of the mutant pool. As 
expected given the strong stress of 100 mM nitrate, this pool was 
delayed by about 92 h in reaching an ODsoo of 0.4 compared with 
only about 40 h in the defined MOLS4 medium without nitrate. 
However, the nitrate pool was harvested and the fitness values 
were determined as previously described (Pels et al., 2013). The 
total number of cells in the final 500 mL culture (1 x 10^' cells) 
was determined by plating for individual colony- forming units, as 
previously described (Pels et al, 2013). The number of doublings 
of the culture was estimated by assuming that only the genes 
with log2 fitness scores >0 (38 genes) contributed significantly 
to the final population. Therefore, the number of unique inser- 
tions in these genes (1904) was considered to be number of cells 
in the starting pool. The complete data from these experiments 
are available at http://desulfovibriomaps.biochem.missouri.edu/ 
fitness/. 

PLASMID AND STRAIN CONSTRUCTION 

Genomic DNA from DvH was isolated with the Wizard® 
Genomic DNA Purification Kit (Promega, USA). Plasmids were 
isolated from both E. coli and DvH with the GeneJET Plasmid 



Miniprep Kit (Fermentas, Thermo Scientific, Glen Burnie, MD). 
All primers were obtained from Integrated DNA Technologies 
(Coralville, lA). The pM03311 and pM03313 plasmids were 
constructed by Sequence- and Ligation-Independent Cloning 
(SLIC) (Li and Elledge, 2007). PGR products from template plas- 
mids were agarose gel-purified to reduce transformation of the 
parent plasmid. All products were cleaned with a Wizard® SV 
Gel and PGR Glean-up System (Promega, USA) before the SLIG 
procedure. The plasmids were constructed by amplification of 
DNA regions (Table S2) with Herculase II polymerase (Agilent 
cat# 600675), as previously described for a similar procedure 
(Parks et al, 2013). DNA products were transformed into Silver 
Efficiency a-Select E. coli cells (Bioline) and plated on solidified 
LG medium (Zane et al., 2010). Electroporation procedures were 
similar to those previously described (Keller et al., 2011) with 
electroporation parameters 1500 V, 250 Q, and 25 |xF. Gells recov- 
ered overnight after electroporation were plated on MOYLS4 
with 1.2 mM thioglycolate as reductant and about 0.2% (wt/vol) 
yeast extract. Sequence confirmation of the mutagenic cassette 
and the complementing gene was performed at the University of 
Missouri DNA Gore Facility (http://www.biotech.missouri.edu/ 
dnacore/). 

NITRATE DETERMINATION 

A scaled-down version of a previously described colorimetric 
method (Gataldo et al, 1975) was used to determine nitrate con- 
centrations. Briefly, 200 [lL of salicylic acid solution (1 g salicylic 
acid dissolved in 20 mL of approximately 98% [vol/vol] sulfuric 
acid) was added to each 25 |xL sample that had been diluted 25- 
fold in deionized water. This was mixed and incubated 20 min at 
room temperature and then 4.75 mL of 2 M NaOH was added to 
each tube. Absorbance at 410 nm was measured with a Genesys 
20 spectrophotometer (Thermo Scientific, Waltham, MA). 
for a standard curve was >0.99, instrument detection limit 
0.1 ±0.1 mM. 

PROTEIN DETERMINATIONS 

Whole cell protein concentrations were determined with the 
Bradford assay (Bradford, 1976) with bovine serum albumin as 
the standard. Absorbance at 595 nm was measured with a Genesys 
20 spectrophotometer (Thermo Scientific, Waltham, MA). 

RESULTS 

RESPONSE OF DvH TO NITRATE EXPOSURE 

It has been reported (Ehas et al, 2009; He et al., 2010a) that 
DvH cells can grow rapidly and abundantly after a long lag phase 
in high nitrate concentrations. It was unclear, however, whether 
this rapid growth was due to elimination of the toxic nitrate, 
some modification of cell metabolism allowing adaptation to the 
continued presence of nitrate, or outgrowth of preexisting nitrate- 
resistant mutants. Furthermore, it was not known whether the 
cells that grew in nitrate had a growth advantage over naive cells 
when subcultured into a fresh medium amended with nitrate. To 
test this, JW710 (Table 1), the parent for making marker exchange 
and markerless deletion strains (Keller et al, 2009), was used. 
JW710 will therefore be referred to as the wild- type control for 
all DvH growth kinetics in this study. This wild-type control was 
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grown in lactate-sulfate medium amended with 100 mM nitrate 
(Figure lA). 

It was determined that, at the end of growth in the presence 
of 100 mM nitrate, no gross consumption of nitrate was detected 
(Table 2). The persistence of the nitrate suggested that the abil- 
ity of DvH to grow in the presence of 100 mM nitrate was due 
to adaptation to nitrate or outgrowth of spontaneous mutants, 
rather than a detoxification of nitrate itself To further confirm 
this lack of nitrate metabolism, nitrite was measured (American 
Public Health Association, 1992) during the lag/inhibition phase 
for JW710 cells exposed to 100 mM nitrate and the nitrite con- 
centration was below 15 ± 5 [lM (data not shown). This nitrite 
concentration is below the concentration (40 (jlM) reported to 
inhibit plated single colonies of D. vulgaris (Haveman et al., 2004) . 
Thus, secondary production of nitrite is not likely the cause of 
nitrate sensitivity. To begin to test whether the nitrate adaptation 
was due to a reversible gene regulation or to successful growth 
of spontaneous mutants in the culture, nitrate-adapted strains 
were subcultured back into medium with or without nitrate. We 
found that nitrate-stressed cultures grew without a prolonged lag 
and maintained nitrate resistance when subcultured into fresh 
medium with 100 mM nitrate (Figure IB). This resistance con- 
tinued over the course of three subcultures (data not shown). 
Further, even nitrate-resistant cultures that were subcultured into 
medium lacking nitrate retained nitrate resistance when subcul- 
tured back into 100 mM nitrate (data not shown). As with the 
original exposure to nitrate, no gross consumption of nitrate was 
observed over the course of the subcultures (Table 2). We sug- 
gest that spontaneous mutations in the culture lead to increased 
nitrate resistance of some cells which then predominate in the 
population. 

FITNESS PROFILING WITH G20 

To test what mutations might be causing this nitrate resistance, 
we employed transposon mutant fitness profiling. A catalogued 
transposon mutant library, which enables high-throughput phe- 
notypic screening, had been generated in G20 prior to that 
produced in DvH (Price et al., 2013). Each mutant strain in 
the library is identified by two unique DNA barcode sequences 
or "tags," the "up" tag and the "down" tag (Oh et al, 2010; 
Deutschbauer et al, 2011). Strain abundance is measured by 
reading the abundance of the barcodes through fluorescence in 




200 



FIGURE 1 I Growth and subculture of wild-type D. vulgaris 
Hildenborough (JW710) in lactate-sulfate medium with inhibitory 
nitrogen species. (A) Growth of DvH with no additions (o) or 100 mM 
sodium nitrate (•) (B) Growth of DvH subcultured from original 100 mM 



microarrays made to detect the barcodes (Pierce et al, 2006, 
2007). 

We predicted that comparison of fitness profiles of nitrite- and 
nitrate-stressed G20 cells would reveal mutants with responses 
unique to nitrate. That is, fitness profiling in the two condi- 
tions would allow us to see which mutants differentially increased 
in relative abundance during a pooled grovrth competition and 
which decreased. The fitness of a particular strain is calculated 
as log2 of the ratio of the relative abundance of the strain after 
growth competition to the relative abundance of the strain before 
growth competition. Therefore, if the relative abundance of a par- 
ticular strain in the pool remained the same before and after 
stress, its fitness would be equal to zero (Oh et al., 2010): 

Fitness of a G20 mutant 

/ barcode microarray signal from cells after stress \ 

= log-, I I 

\ barcode microarray signal from cells prior to stress / 



Table 2 | Nitrate concentrations from stationary phase cultures of 
D. vulgaris Hildenborough grown in M0LS4 medium amended as 
indicated. 



Subculture^ 


Inoculum history 


Amendment 


[NOj]'' in mM 




Lactate/SOj" 


none 


not measured 


0-B'= 


Lactate/SO^" 


lOOmM NO" 


101 ±3 


l-A-^ 


From 0-B 


none 


8 ± 1 


1-6=^ 


From 0-B 


lOOmM NO" 


102 ± 3 


2-A 


From 1-B 


none 


6±2 


2-B 


From 1-B 


lOOmM NO3 


99 ± 1 


3-A 


From 2-A 


none 


not detected 


3-B 


From 2-A 


lOOmM NO3 


101 ±7 


3-C 


From 2-B 


none 


8±2 


3-D 


From 2-B 


lOOmM NO3 


103 ± 5 



"Inocula were 6.5% (vol/vol). 

''Concentrations determined from triplicate determinations with standard devia- 
tions shown. 

Growth curves in Figure 1A. 

Growth curves in Figure IB. 




nitrate culture shown in (A) no additions (o), 100 mM sodium nitrate (•). 
Approximately 6-6.5% (vol/vol) inocula were used for the original culture and 
subcultures. Readings reflect averages of three samples, and errors bars 
show standard deviations. 
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If a strain decreased in relative abundance after the stress con- 
dition because it was outcompeted or unable to cope with the 
stress, it would have negative ( <0) fitness. If its relative abundance 
increased, it would have positive (>0) fitness. For the pools, the 
fitness calculated for a particular gene, referred to as the "mean 
log ratio," is expressed as log2 of the average fitness of strains with 
a mutation in that particular gene. 

The G20 pools were grown in lactate -sulfate medium amended 
with 150 mM NaNOs, 150 mM KNO3, 150 mM NaCl, 150 mM 
KCl, or 0.25 mM NaN02. NaCl and KCl conditions were osmotic 
controls; KNO3 vs. KCl allowed a control for anion specificity. 
Concentrations of 150 mM nitrate and 0.25 mM nitrite were cho- 
sen because these concentrations severely but not completely 
inhibited wild-type G20 (Figures 2 A,B). Since a long lag phase 
had been observed before exponential growth of nitrate-stressed 



cultures of both DvH and G20 in lactate-sulfate conditions 
(Figures lA, 2A), it was reasonable to hypothesize that sponta- 
neous mutants in wild-type cultures were selected in the pop- 
ulation after the lag. Therefore, this hypothesis was confirmed 
when we found that, in the transposon mutant pools, several 
mutant strains predominated in cultures growing in the pres- 
ence of nitrate (Table 3; Table S3). That is, transposon insertion 
conferred a growth advantage, and therefore a high fitness, on 
these particular strains in the nitrate stress condition. The top 10 
genes interrupted in the strains that grew abundantly in sodium 
nitrate had fitness values (mean log ratios) greater than 2 in 
that condition, but fitness values less than 0.25 in sodium chlo- 
ride, potassium chloride, sodium nitrite, and in the absence of 
stress (Table S3). In contrast, interruption of those same genes 
was advantageous in both sodium nitrate and potassium nitrate 




50 75 
Time (h) 



FIGURE 2 I Growth of wild-type D. alaskensis G20 in lactate-sulfate 
medium with inhibitory nitrogen species. (A) Growth of G20 with no 
additions (o), 150 mM sodium nitrate (□), or 150 potassium nitrate (■). (B) 



125 



Growth of G20 with no additions (o) or 0.25 mM sodium nitrite (•). 
Approximately 4.5% (vol/vol) inocula were used. Readings reflect averages of 
three samples, and errors bars show standard deviations. 



Table 3 | Desulfovibrio genes interrupted in strains with high fitness in lactate-sulfate conditions amended with sodium nitrate. 
G20 gene (Dde) fitness^ DvH homolog (DVU) fitness'' Annotations 



2702 

0597 

0598 

0600 
0601 

0602 

0603 
0604 

0605 



4.23 

2.25 

3.01 

2.88 
3.51 



0916 

no homolog 

0251 

0250 
0249 



2.34 0248 (pseudo-gene) 

3.08 0247 

3.41 0246 

2.08 0245 



3.81 



No data 



11.44 



-5.82 
3.86 



1.43 



9.14 



2.18 



-6.04 



ATrich DNA binding protein (COG2344); Transcriptional repressor, redox-sensing, 
Rex (IPR022876) 

Uncharacterized protein conserved in archaea (COG2043); Protein of unknown 
function DUF169 (IPR003748) 

Transmembrane protein TauE like (IPR002781); predicted permease (COG0730 ); 
sulfite exporter TauE/SafE (pfam01925) 
Conserved hypothetical protein 

PtxB, putative (http://microbesonline.org/); ABC-type phosphate/phosphonate 
transport system, periplasmic component (COG3221); outer 
membrane-associated homodimer (Walian et a!., 2012) 
Signal transduction histidine kinase (COG5002); PAS fold (IPR013767); 
ATPase-like, ATP-binding domain (IPR003594); HAMP linker domain (IPR003660); 
PAC motif (IPR001610) 

Signal transduction response regulator, receiver domain (IPR001789); CheY-like 

superfamily (IPR011006); ntrX (http://microbesonline.org/) 

Pyruvate phosphate dikinase, PEP/pyruvate-binding (IPR002192); PEP-utilizing 

enzyme, mobile domain (IPR008279); ATP-grasp fold, subdomain 1 (IPR013815); 

ATP-grasp fold, subdomain 2 (IPR013816) 

Protein-tyrosine/Dual-specificity phosphatase (IPR000387) 



I092 ( £:^^;:7J^[:r.TJ"^. ) : ff^^^ of stationary-phase G20 cultures grown for about 3.3 doublings (about 63h) in lactate-sulfate medium 



^ barcode microanay signal from cells prior to stress 
amended with 150 mM sodium nitrate. 

b 



''log2[[* 



insertions in c 



length of genT" ) / ( * 'length "of aii%nlT'' ) ] f'f"ess determined from mid-log phase DvH cultures grown for about 25.5 doublings (92 h) in lactate-sulfate 



medium amended with 100 mM nitrate. 
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(Table S3). Growth of the mutants in both sahs of nitrate sup- 
ports the specificity of the nitrate anion as the driver for selection 
of these resistant mutants. The interrupted gene conferring the 
highest fitness in sodium nitrate was Dde_2702 (Table 3), a gene 
annotated as encoding Rex, a redox-sensing regulatory protein 
(Ravcheev et al, 2012). Particularly surprising was the high fit- 
ness conferred by mutation of a cluster of poorly annotated genes, 
Dde_0597 through Dde_0605, hereafter called the "nitrate clus- 
ter." Both the rex gene and the nitrate cluster (Table 3) have 
homologs in DvH. Because of these homologies, including shared 
synteny of the nitrate cluster in G20 (Figure 3), it seemed reason- 
able that mutations of the homologs in DvH would confer similar 
nitrate-resistant phenotypes. 

FITNESS PROFILING WITH DvH 

In order to test the hypothesis that both G20 and DvH used 
the same mechanisms for nitrate resistance, we had the oppor- 
tunity to employ a different high-throughput fitness profiling 
method, Transposon Liquid Enrichment sequencing (TnLE-seq). 
This method (Pels et al., 2013) is based on deep sequencing of 
random transposon mutations to query DvH. TnLE-seq is a mod- 
ification of the HITS (High-throughput Insertion Tracking by 
deep Sequencing) (Gawronski et al., 2009), Tn-seq (Van Opijnen 
et al., 2009), and TraDIS (Transposon Directed Insertion-site 
Sequencing) (Langridge et al, 2009) methods. However, the 
TnLE-seq method developed for DvH is especially well-adapted 
to oxygen-sensitive bacteria that have low electroporation effi- 
ciency (Pels et al., 2013). The mutated culture is grown in control 
vs. stress conditions, and deep sequencing then determines the 
abundance and location of mutations at the end of growth. 
Because of the differences in methods, the calculation of fitness 
is also different from that of the mutant library experiment (Pels 
et al., 2013). The fitness value shown below is in logaR format, for 
easier comparison with the G20 pools: 



Pitness of a DvH gene = log2 



/ 



# insertions in gene \ 
length of gene / 

# insertions in all genes \ 
length of all genes / 



As previously described (Pels et al., 2013), fitness was calculated 
from insertions only in the 5-85% region of the coding sequence 
of genes, as such insertions are more likely to impair the func- 
tion of gene products. As with the G20 pool described above, 
negative fitness indicates a fitness defect, whereas positive fitness 
indicates that the mutation confers a fitness advantage in that par- 
ticular condition. Por nitrate stress, the transposon mutants were 



grown in lactate-sulfate medium amended with 100 mM sodium 
nitrate. As expected, the results showed that mutations in a pre- 
dicted rex gene annotated as encoding a transcriptional regulator 
(DVU0916) as well as mutations in homologs of the G20 "nitrate 
cluster" (DVU0251, DVU0249, DVU0247, and DVU0246) con- 
ferred fitness values among the ten highest values (Table 3; Table 
S4). In fact, mutation of DVU0251 led to the highest fitness value, 
11.44, or 2780-fold. Essentially, there was a "jackpot effect" in 
which a small percentage of mutants predominated in the pop- 
ulation, a consistent result between the DvH and G20 fitness 
experiments. Despite these consistencies, in-depth, individual 
mutant analysis was necessary to confirm and elucidate the results 
of high-throughput fitness profiling (Deutschbauer et al, 2011). 

CONFIRMATION OF FITNESS PROFILING WITH INDIVIDUAL MUTANTS 

DvH was chosen for confirmation studies because a catalogued 
transposon mutant library of DvH was also available, in-frame 
deletion mutants can be made with greater facility (Keller et al., 
2009), and complementation of mutants is readily accomplished. 
For an initial confirmation of the physiological relevance of 
the "nitrate cluster" to nitrate resistance, we determined growth 
kinetics of five DvH isolated mutants with transposon inser- 
tions in genes in that cluster (Figures 4A-E). The control strain 
used had a transposon at an intergenic position 327 base pairs 
upstream of a gene encoding a "random" hypothetical pro- 
tein, DVU0590, and, therefore, was not predicted to be involved 
in nitrate stress responses. We found that the mutants with 
transposon insertions in DVU0246, DVU0247, DVU0250, and 
DVU0251 grew with indistinguishable kinetics with or with- 
out 100 mM nitrate. In contrast, 100 mM nitrate inhibited the 
control and the DVU0245 transposon mutants (Figures 4 A-F). 
Inhibition of the DVU0245 mutant in high nitrate is consis- 
tent with the low fitness of the DVU0245 mutant in nitrate 
(fitness —6.04, Table 3) and consistent with the inhibition of a 
deletion mutant of DVU0245 (data not shown). None of the 
mutants grew better than the control strain in 1 mM sodium 
nitrite (Figures 4A-F). We interpret these results to mean that 
the growth advantage of these mutants is specific to nitrate and 
not simply an advantage in the presence of inhibitory nitrogen 
species. 

The results of the G20 and DvH fitness studies indicated sim- 
ilar responses of homologous genes. In the G20 results, the pre- 
dicted rex mutant had the highest fitness (Dde_2702). The fitness 
score for the DvH rex mutant was also in the top ten fitness scores, 
along with four DvH "nitrate cluster" (DVU0251, DVU0247, 
DVU0249, DVU0246) mutants (Table 3; Table S2). However, 
both the DVU0245 and the DVU0250 mutants had low fitness 
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FIGURE 3 1 Desulfovibrio nitrate resistance gene cluster. Operon predictions were from http://microbesonline.org/; boxes represent predicted genes, arrows 
indicate direction of transcription, and contiguous boxes ending in an arrow represent predicted operons. 
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FIGURE 4 I Growth of D. vulgaris Hildenborough transposon mutants presence of no additions (o), 100 mM sodium nitrate (□), or 1 mM 

in lactate-sulfate medium with inhibitory nitrogen species. Growth sodium nitrite (•). Approximately 5.3% (vol/vol) inocula were used, 

of the DVU0245 (A), DVU0246 (B), DVU0247 (C), DVU0250 (D), Optical density readings reflect averages of three samples, and errors 

DVU0251 (E), and intergenic transposon control (F) mutants in the bars show standard deviations. 



in nitrate, whereas mutants of their G20 homologs had high fit- 
ness in nitrate (Table 3). The nitrate resistance of the pooled vs. 
the isolated DvH mutants was also not entirely consistent. The 
rapid growth of the isolated DVU0250 mutant (Figure 4D) in the 
presence of nitrate was unexpected because of the low TnLE-seq 
fitness conferred in the same condition by mutation of DVU0250 
(Tables). In addition, preliminary growth kinetic data (not 
shown) suggest that the transposon mutant of DVU0248, which 
is annotated as a pseudogene in DvH, has little or no growth 
advantage over the control strain in the presence of 100 mM 
nitrate. In contrast, in the pooled experiment the DVU0248 
mutations conferred positive fitness in nitrate (Table 3). While 
more data wiU be needed either to confirm or to change the 
annotation of DVU0248 as a pseudogene, we suggest that simi- 
lar nitrate-resistance mechanisms are operating in G20 and DvH. 
The discrepancies between pooled and individual mutant stud- 
ies confirm the need for follow-up studies of high-throughput 
experiments. 



Such follow-up was pursued with gene deletion and comple- 
mentation of the DvH gene encoding the predicted transcription 
regulator Rex. Interruption of the rex gene conferred the highest 
fitness in G20 but not in DvH (Table 3). We found that a deletion 
of DvH rex (DVU0916) had a clear advantage over the JW710 par- 
ent strain in lactate-sulfate medium with 100 mM added nitrate 
(Figures 5 A,B)- Like the "nitrate cluster" transposon mutants 
described above (Figures 4A-F), the rex mutant is not demon- 
strably resistant to nitrite (Figures 5A,B). Interestingly, when the 
mutant was complemented with rex expressed from a constitutive 
promoter, the phenotype in the presence of 100 mM nitrate was 
different from either parent or mutant phenotypes (Figure 5C). 
In contrast, the parent strain with rex overexpressed appeared to 
be at least as sensitive to nitrate as the parent strain (Figure 5C). 
The unique phenotype of the complemented mutant may result 
from some of the bacterial population losing the plasmid contain- 
ing the complemented rex gene, in spite of antibiotic selection. 
While spectinomycin selects for plasmid retention, nitrate should 
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50 Time (h) lOO 



FIGURE 5 I Growth of DvH wild-type vs. Arex mutant in lactate-sulfate 
medium with inhibitory nitrogen species. (A,B) Growth of the "wild-type " 
parental strain JW710 (A) vs. Arex (B) mutant in the presence of no 
additions (o), 100 mM sodium nitrate (□), or 1 mM sodium nitrite (•). 
Approximately 6% (vol/vol) inocula were used. (C,D) show four strains grown 
with 100 mM nitrate (C) or no additions (D). Wild-type with empty vector 



50 Time{h) 100 



[JW710(pMO9075)l,(A); wild-type with rex overexpression plasmid 
[JW710(pMO3313)l, (A); Arex strain with empty vector lJW3311(pMO9075)l, 
(x); Arex strain with rex complement plasmid 1JW3311 (pM03313)l, (■). 
Approximately 75% (vol/vol) inocula were used. Readings reflect averages of 
three samples, and errors bars, which were often within the symbols, show 
standard deviations. 



select for plasmid loss in a Arex strain grown in high nitrate. Cells 
containing the plasmid may produce enough of the antibiotic- 
modifying enzyme to confer sufficient resistance to allow other 
cells to survive without containing the plasmid. If this is the 
case, then the Arex cells containing the plasmid should grow 
slowly while those which have lost the plasmid should grow more 
rapidly in the presence of 100 mM nitrate. The result would 
be a population growth rate in-between that seen for wild-type 
vs. Arex strains. Indeed, the phenotype of the complemented 
Arex strain exhibits this growth property (Figure 5C). Finally, 
nitrate concentrations in the cultures with empty vector or rex 
complementing plasmids (Figures 5C,D) were measured colori- 
metrically at the end of growth. As with the wild-type cultures 
described above, gross consumption of nitrate was not detected 
for any of these strains (data not shown). This is evidence of 
genuine nitrate resistance in these strains. Taken together, these 
growth and gene fitness data support transcriptomic predictions 
that nitrate stress responses involve mechanisms independent of 
nitrite stress responses (He et al, 2010a). 

DISCUSSION 

Whereas Elias et al. (2009) and Bender et al. (2007) reported 
mutations that led to increased sensitivity to both nitrate and 
nitrite, here we report the unexpected discovery of DvH mutants 
with increased resistance to nitrate but not nitrite. The data pre- 
sented confirm that the rex deletion and the "nitrate cluster" 
transposon mutants, the top candidates from fitness profiling, 
confer resistance to nitrate in DvH. Such resistance also devel- 
oped in the non-mutagenized DvH parental strain after subcul- 
ture from 100 mM nitrate (Figure IB), likely as a result of the 



outgrowth of preexisting spontaneous nitrate-resistant mutants. 
The lack of nitrate metabolism of DvH is consistent with a report 
that 10 mM nitrate did not noticeably inhibit DvH (Haveman 
et al., 2004). Furthermore, DvH has been shown to reduce nitrite 
(Haveman et al., 2004), but not nitrate. 

The ability of mutations in a subset of non-essential genes 
to confer nitrate resistance may in part account for the recently 
reported fluctuating sulfide levels produced by SRB in a biore- 
actor inoculated with oU from a Canadian oil field (Callbeck 
et al., 2013). In this bioreactor, sulfide production was completely 
inhibited during pulses of 100 mM nitrate. However, after each 
pulse, sulfide production resumed, indicating that some SRB per- 
sisted in the presence of the nitrate (Callbeck et al, 2013). The 
results of the work presented here suggest that persistence of 
SRB in nitrate-treated oil reservoirs may be the result of mutant 
resistance. Even if total oil-well nitrate concentrations reach low 
millimolar levels, the initial concentration of nitrate near the 
injection site will be much higher than this. For example, the peak 
nitrate concentration in one study of pulsed nitrate injection was 
reported as 760 mM (Voordouw et al, 2009). Resistance to nitrate 
in the presence of a mixed culture is consistent with preliminary 
fitness profiling data from G20 grown in coculture with the nitrate 
reducer Pseudomonas stutzeri RCH2 in the presence of 100 mM 
nitrate. Under these mixed-culture conditions, the "nitrate clus- 
ter" mutant and the rex mutant gained a fitness advantage (A. 
Deutschbauer, unpublished data) very similar to that observed in 
pools of G20 mutants alone in the presence of 150 mM nitrate 
(Tables). 

These fitness studies bring clarity to questions that neither 
transcriptomic nor proteomic data could answer. While "omics" 
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Studies can assist detection and monitoring of changes in the 
metabolism of bacteria in contaminated environments (Steinberg 
et al., 2008), they are not sufficient (Torres-Garcla et al, 2009) 
for elucidating underlying inhibitory mechanisms. In fact, there 
are poor correlations between the expression of transcripts and 
the expression of proteins in DvH in response to nitrate stress 
(Redding et al, 2006; He et al, 2010a). He et al. (2010a) reported 
28 genes for which the mRNA and protein levels were both sig- 
nificantly changed in nitrate stress conditions. However, for 7 
of these 28 genes, the mRNA was significantly downregulated 
while the protein was significantly upregulated (He et al., 2010a). 
Although this poor correlation may be a result of meaningful 
regulatory mechanisms (Lu et al., 2007), transcript abundance 
is difficult to interpret and does not always correlate well with 
gene fitness (Price et al., 2013). For example, there is an upward 
trend of expression of the "nitrate cluster" genes in both nitrate 
(He et al., 2010a) and nitrite (He et al, 2006) stress conditions 
(http://microbesonline.org/). Because mutation of the nitrate 
cluster genes confers a growth advantage in high nitrate, increased 
expression of these genes should be detrimental to growth of 
DvH in high nitrate. This is consistent with the recent deduction 
(Price et al., 2013) that, counterintuitively, detrimental bacterial 
genes are often not downregulated. This apparent suboptimal 
regulation of the nitrate cluster genes in the presence of high con- 
centrations of nitrate likely contributes to the explanation of why 
their interruption confers such a strong growth advantage. 

The roles of the genes of the nitrate cluster in nitrate sensi- 
tivity are not immediately obvious from their poor annotation 
(Tables). The native functions of the nitrate cluster genes are 
not likely involved with nitrate, since neither DvH (Seitz and 
Cypionka, 1986; Pereira et al, 2000) nor G20 (Wall, unpubhshed 
data) have been shown to use nitrate for energy conservation. We 
hypothesize that they may allow non-specific nitrate transport by 
a leaky thiosulfate transporter. The G20 mutants in this cluster 
were mildly sick when grown with 10 mM thiosulfate as a termi- 
nal electron acceptor (A. Deutschbauer, unpublished data), which 
suggests a possible role for these genes in thiosulfate uptake. 
Preliminary data also show that the Arex DvH strain described 
here grows more slowly than the parent strain with 30 mM thio- 
sulfate as a terminal electron acceptor (Christensen, unpublished 
data). Therefore, we suggest that mutation of these genes might 
relieve nitrate inhibition by barring entry of nitrate into the cell. 
Further study will be needed to explore their native functions. 

Several additional genes outside of the "nitrate cluster" appear 
to have high fitness values during nitrate stress, indicating that 
their absence may improve growth. Follow-up studies with indi- 
vidual mutants will be necessary to confirm these predictions. The 
results from this study clearly indicate that DvH and G20 have 
common nitrate resistance mechanisms that should be considered 
in environmental modeling. 

AUTHOR CONTRIBUTIONS 

Hannah L. Korte, Samuel R. Fels and Judy D. Wall designed 
the experiments. Samuel R. Fels conducted the TnLE-seq exper- 
iments and analyzed them. Geoff A. Christensen constructed the 
Arex mutant and showed its impairment on thiosulfate. Morgan 
N. Price and Adam M. Deutschbauer generated and analyzed 



the G20 fitness profiling data. Jennifer V. Kuehl constructed 
the G20 transposon mutant library. Grant M. Zane constructed 
the DvH transposon mutant library, confirmed and provided 
the individual transposon mutants tested. Adam P. Arkin and 
Judy D. Wall supervised the project. Hannah L. Korte, Morgan 
N. Price, and Judy D. Wall interpreted the data and wrote the 
manuscript. 

ACKNOWLEDGMENTS 

This work, conducted by ENIGMA- Ecosystems and Networks 
Integrated with Genes and Molecular Assemblies (http://enigma. 
lbl.gov/), a Scientific Focus Area Program at Lawrence Berkeley 
National Laboratory, was supported by the Office of Science, 
Office of Biological and Environmental Research, of the 
U.S. Department of Energy under Contract No. DE-AC02- 
05CH11231. 

SUPPLEMENTARY MATERIAL 

The Supplementary Material for this article can be found online 
at: http://www.frontiersin.org/)ournal/10.3389/fmicb.2014. 
00153/abstract 

REFERENCES 

American Public Health Association. (1992). Standard Methods for The 
Examination of Water and Waste Water, 18th Edn. Washington, DC: American 
Water Works Association and Water Pollution Control Federation. 

Bender, K. S., Yen, H. C. B., Hemme, C. L., Yang, Z., He, Z., He, Q., 
et al. (2007). Analysis of a ferric uptake regulator (Fur) mutant of 
Desulfovibrio vulgaris Hildenborough. Appl. Environ. Microbiol. 73, 5389-5400. 
doi: 10.1128/AEM.00276-07 

Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of 
microgram quantities of protein utilizing the principle of protein-dye binding. 
Anal. Biochem. 72, 248-254. doi: 10.1016/0003-2697(76)90527-3 

Caffrey, S. M., and Voordouw, G. (2010). Effect of sulfide on growth physiol- 
ogy and gene expression of Desulfovibrio vulgaris Hildenborough. Antonie Van 
Leeuwenhoek 97, 11-20. doi: 10.1007/sl0482-009-9383-y 

Callbeck, M. C, Agrawal, A., and Voordouw, G. (2013). Acetate production from 
oil under sulfate-reducing conditions in bioreactors injected with sulfate and 
nitrate. App/. Environ. Microbiol. 79, 5059-5068. doi: 10.1128/AEM.01251-13 

Cataldo, D. A., Haroon, M., Scharder, L. E., and Youngs, V. L. (1975). Rapid 
colorimetric determination of nitrate in plant tissue by nitration of sali- 
cylic acid. Comm. Soil Sci. Plant Anal. 6, 71-80. doi: 10.1080/0010362750 
9366547 

Chardin, B., Dolla, A., Chaspoul, R, Fardeau, M., Gallice, P., and Bruschi, M. 
(2003). Bioremediation of chromate: thermodynamic analysis of the effects of 
Cr(VI) on sulfate-reducing bacteria. Appl. Microbiol. Biotech. 60, 352-360. doi: 
10.1007/S00253-002-1091-8 

Chhabra, S. R., He, Q., Huang, K. H., Gaucher, S. P., Aim, E. J., He, Z., et al. (2006). 
Global analysis of heat shock response in Desulfovibrio vulgaris Hildenborough. 
/. Bacteriol. 188, 1817-1828. doi: 10.1128/JB.188. 5. 1817-1828.2006 

Clark, M. E., He, Q., He, Z., Huang, K. H., Aim, E. J., Wan, X. R, et al. (2006). 
Temporal transcriptomic analysis as Desulfovibrio vulgaris Hildenborough tran- 
sitions into stationary phase during electron donor depletion. Appl. Environ. 
Microbiol. 72, 5578-5588. doi: 10.1128/AEM.00284-06 

Deutschbauer, A., Price, M. N., Wetmore, K. JVI., Shao, W, Baumohl, J. K., Xu, Z., 
et al. (2011). Evidence-based annotation of gene function in Shewanella onei- 
densis MR-1 using genome-wide fitness profiling across 121 conditions. PLoS 
Genet. 7:el002385. doi: 10.1371/journal.pgen.l002385 

EUas, D. A., Mukhopadhyay, A., Joachimiak, M. P., Drury, E. C., Redding, A. 
M., Yen, H. C. B., et al. (2009). Expression profiling of hypothetical genes in 
Desulfovibrio vulgaris leads to improved functional annotation. Nucleic Acids 
Res. 37, 2926-2939. doi: 10.1093/nar/gkpl64 

Fels, S. R., Zane, G. M., Blake, S. M., and Wall, J. D. (2013). Rapid Transposon 
Liquid Enrichment Sequencing (TnLE-seq) for gene fitness evaluation in 



Frontiers in Microbiology | Microbial Physiology and Metabolism 



April 2014 | Volume 5 | Article 153 | 10 



Korte et al. 



Desulfovibrio resistance to high nitrate 



underdeveloped bacterial systems. Appl Environ. Microbiol 79, 7510-7517. doi: 
10.1128/AEM.02051-13 

Gawronski, J. D., Wong, S. M. S., Giannoukos, G., Ward, D. V., and Akerley, B. 
J. (2009). Tracking insertion mutants within libraries by deep sequencing and 
a genome-wide screen for Haemophilus genes required in the lung. Proc. Natl 
Acad. Set. U.S.A. 106, 16422-16427. doi: 10.1073/pnas.0906627106 

Green, S. J., Prakash, O., Jasrotia, P., Overholt, W. A., Cardenas, E., Hubbard, D., 
et al. (2012). Denitrifying bacteria from the genus rhodanobacter dominate bac- 
terial communities in the highly contaminated subsurface of a nuclear legacy 
waste site. App/. Environ. Microbiol 78, 1039-1047. doi: 10.1128/AEM.06435-11 

Grigoryan, A. A., Cornish, S. L., Buziak, B., Lin, S., Cavallaro, A., Arensdorf J. J., 
et al. (2008). Competitive oxidation of volatile fatty acids by sulfate- and nitrate- 
reducing bacteria from an oil field in Argentina. App/. Environ. Microbiol 74, 
4324-4335. doi: 10.1128/AEM.00419-08 

Hauser, L. J., Land, M. L., Brown, S. D., Larimer, F., Keller, K. L., Rapp- 
Giles, B. J., et al. (2011). Complete genome sequence and updated anno- 
tation of Desulfovibrio alaskensis G20. /. Bacteriol 193, 4268-4269. doi: 
10.1128/JB.05400-11 

Haveman, S. A., Greene, E. A., Stilwell, C. P., Voordouw, J. K., and Voordouw, 
G. (2004). Physiological and gene expression analysis of inhibition of 
Desulfovibrio vulgaris HUdenborough by nitrite. /. Bacteriol 186, 7944—7950. 
doi: 10.1128/JB.186.23.7944-7950.2004 

Haveman, S. A., Greene, E. A., and Voordouw, G. (2005). Gene expression analy- 
sis of the mechanism of inhibition of Desulfovibrio vulgaris HUdenborough by 
nitrate-reducing, sulfide-oxidizing bacteria. Environ. Microbiol 7, 1461-1465. 
doi: 10.1 1 1 l/j.l462-2920.2005.00834.x 

Hazen, T. C, and .Stahl, D. A. (2006). Using the stress response to monitor process 
control: pathways to more effective bioremediation. Curr. Opin. Biotechnol. 17, 
285-290. doi: 10.1016/j.copbio.2006. 03.004 

He, Q., He, Z., Joyner, D. C, Joachimiak, M., Price, M. N., Yang, Z. K., et al. {2010a). 
Impact of elevated nitrate on sulfate-reducing bacteria: a comparative study of 
Desulfovibrio vulgaris. ISME J. 4, 1386-1397. doi: 10.1038/ismej.2010.59 

He, Q., Huang, K. H., He, Z., Aim, E. J., Fields, M. W, Hazen, T. C, et al. (2006). 
Energetic consequences of nitrite stress in Desulfovibrio vulgaris HUdenborough, 
inferred from global transcriptional analysis. Appl Environ. Microbiol 72, 
4370^381. doi: 10.1128/AEM.02609-05 

He, Z., Zhou, A., Baidoo, E., He, Q., Joachimiak, M. R, Benke, P., et al. (2010b). 
Global transcriptional, physiological and metabolite analyses of the responses 
of Desulfovibrio vulgaris Hildenborough to salt adaptation. Appl Environ. 
Microbiol 76, 1574-1586. doi: 10.1128/AEM.02141-09 

Heidelberg, J. P., Seshadri, R., Haveman, S. A., Hemme, C. L., Paulsen, I. T., 
Kolonay, J. P., et al. (2004). The genome sequence of the anaerobic, sulfate- 
reducing bacterium Desulfovibrio vulgaris Hildenborough. Nat Biotech. 22, 
554-559. doi: 10.1038/nbt959 

Jalali, K., and Baldwin, S. A. (2000). The role of sulphate reducing bacteria in 
copper removal from aqueous sulphate solutions. Water Res. 34, 797—806. doi: 
10.1016/50043-1354(99)00194-3 

Jiang, W., and Fan, W. (2008). Bioremediation of heavy metal-contaminated 
soils by sulfate-reducing bacteria. Ann. N.Y. Acad. Sci. 1140, 446-454. doi: 
10.1196/annals.l454.050 

Keller, K. L., Bender, K. S., and WaU, J. D. (2009). Development of a markerless 
genetic exchange system for Desulfovibrio vulgaris Hildenborough and its use 
in generating a strain with increased transformation efficiency. Appl Environ. 
Microbiol 75, 7682-7691. doi: 10.1128/AEM.01839-09 

Keller, K. L., Rapp-Giles, B. J., Semkiw, E. S., Porat, I., Brown, S. D., and Wall, J. 
D. (2014). New model for electron flow for sulfate reduction in Desulfovibrio 
alaskensis G20. Appl Environ. Microbiol 80, 855-868. doi: 10.1128/AEM. 
02963-13 

Keller, K. L., WaU, J. D., and Chhabra, S. (2011). Methods for engineermg sulfate 
reducing bacteria of the genus Desulfovibrio. Methods Enzymol. 497, 503-517. 
doi: 10.1016/B978-0-12-385075-1.00022-6 

Kircher, M., Sawyer, S., and Meyer, M. (2012). Double indexing overcomes inac- 
curacies in multiplex sequencing on the lUumina platform. Nucleic Acids Res. 
40:e3. doi: 10.1093/nar/gla-771 

Langridge, G. C., Phan, M. D., Turner, D. J., Perkins, T. T, Parts, L., Haase, J., 
et al. (2009). Simultaneous assay of every Salmonella typhi gene using one 
million transposon mutants. Genome Res. 19, 2308-2316. doi: 10.1101/gr.0970 
97.109 



Larsen, R. A., WUson, M. M., Guss, A. M., and Metcalf, W. W. (2002). Genetic 
analysis of pigment biosynthesis in Xanthobacter autotrophicus Py2 using a 
new, highly efficient transposon mutagenesis system that is functional in a 
wide variety of bacteria. Arch. Microbiol 178, 193-201. doi: 10.1007/s00203-00 

2-0442-2 

Li, M. Z., and EUedge, S. J. (2007). Harnessing homologous recombination 

in vitro to generate recombinant DNA via SLIC. Nat. Methods 4, 251-256. doi: 

10.1038/nmethlOlO 

Ligthelm, D. J., De Boer, R. B., Brint, J. R, and Schulte, W. M. (1991). "Reservoir 
souring. An analytical model for H2S generation and transportation in an 
oil reservoir owing to bacterial activity," in Offshore Europe 91 — Proceedings, 
(Aberdeen, Scotland: Society of Petroleixm Engineers of AIME), 369-378. doi: 
10.2118/23141-MS 

Lloyd, J. R., Ridley, J., Khizniak, T., Lyalikova, N. N., and Macaskie, L. E. (1999). 
Reduction of technetium by Desulfovibrio desulfuricans: biocatalyst charac- 
terization and use in a flowthrough bioreactor. Appl Environ. Microbiol 65, 
2691-2696. 

Lovley, D. R., Roden, E. E., PhUUps, E. J. P., and Woodward, J. C. (1993a). Enzymatic 
iron and uranium reduction by sulfate-reducing bacteria. Marine Ceol. 113, 
41-53. doi: 10.1016/0025-3227(93)90148-0 

Lovley, D. R., Widman, R K., Woodward, J. C., and PhiUips, E. J. R (1993b). 
Reduction of uranium by cytochrome C3 of Desulfovibrio vulgaris. Appl Environ. 
Microbiol 59, 3572-3576. 

Lu, R, Vogel, C, Wang, R., Yao, X., and Marcotte, E. M. (2007). Absolute 
protein expression profiling estimates the relative contributions of transcrip- 
tional and translational regulation. Nat. Biotech. 25, 117-124. doi: 10.1038/ 
nbtl270 

Martins, M., Faleiro, M. L., Barros, R. J., Verissimo, A. R., Barreiros, M. A., and 
Costa, M. C. (2009). Characterization and activity studies of highly heavy metal 
resistant sulphate-reducing bacteria to be used in acid mine drainage decontam- 
ination. /. Hazard. Mater 166, 706-713. doi: 10.1016/j.jhazmat.2008.11.088 

McCready, R. G. L., Gould, W. D., and Cook, F. D. (1983). Respiratory 
nitrate reduction by Desulfovibrio sp. Arch. Microbiol. 135, 182-185. doi: 
10.1007/BF00414476 

Mukhopadhyay A., He, Z., Aim, E. J., Arkin, A. R, Baidoo, E. E., Borglin, S. C, 
et al. (2006). Salt stress in Desulfovibrio vulgaris HUdenborough: an integrated 
genomics approach./. Bacteriol 188,4068^078. doi: 10.1128/JB.01921-05 

Mukhopadhyay, A., Redding, A. M., Joachimiak, M. P., Arkin, A. P., Borglin, 
S. E., Dehal, P. S., et al. (2007). CeU-wide responses to low-oxygen expo- 
sure in Desulfovibrio vulgaris HUdenborough. /. Bacteriol. 189, 5996-6010. doi: 
10.1128/JB.00368-07 

Oh, J., Fung, E., Price, M. N., Dehal, R S., Davis, R. W, Giaever, G., et al. (2010). A 
universal tagmodule collection for paraUel genetic analysis of microorganisms. 
Nucleic Acids Res. 38, 14. doi: 10.1093/nar/gkq419 

Parks, J. M., Jobs, A., Podar, M., Bridou, R., Hurt, R. A. Jr., Smith, S. D., et al. (2013). 
The genetic basis for bacterial mercixry methylation. Science 339, 1332—1335. 
doi: 10.1126/science.l230667 

Pereira, I. A., LeGall, J., Xavier, A. V., and Teixeira, M. (2000). Characterization 
of a heme c nitrite reductase from a non-ammonifying microorganism, 
Desulfovibrio vulgaris Hildenborough. Biochim. Biophys. Acta 1481, 119—130. 
doi: 10.1016/S0167-4838(00)00111-4 

Pierce, S. E., Davis, R. W., Nislow, C, and Giaever, G. (2007). Genome-wide analysis 
of barcoded Saccharomyces cerevisiae gene-deletion mutants in pooled cultures. 
Nat Protoc. 2, 2958-2974. doi: 10.1038/nprot.2007.427 

Pierce, S. E., Fung, E. L., JaramiUo, D. E, Chu, A. M., Davis, R. W, Nislow, C, 
et al. (2006). A unique and universal molecular barcode array. Nat. Methods 3, 
601-603. doi: 10.1038/nmeth905 

Postgate, J. R. (1984). The Sulfate-Reducing Bacteria, 2nd Edn. Cambridge: 
Cambridge University Press. 

Price, M. N., Deutschbauer, A. M., Skerker, J. M., Wetmore, K. M., Ruths, T., 
Mar, J. S., et al. (2013). Indirect and suboptimal control of gene expression is 
widespread in bacteria. Mol Syst. Biol 9, 660. doi: 10.1038/msb.2013.16 

Ravcheev, D. A., Li, X., Latif H., Zengler, K., Leyn, S. A., Korostelev, Y. D., et al. 
(2012). Transcriptional regulation of central carbon and energy metabolism in 
bacteria by redox-responsive repressor rex. /. Bacteriol 194, 1145-1157. doi: 
10.1128/JB.06412-11 

Redding, A. M., Mukhopadhyay, A., Joyner, D. C, Hazen, T. C, and Keasling, 
J. D. (2006). Study of nitrate stress in Desulfovibrio vulgaris HUdenborough 



www.frontiersin.org 



April 2014 I Volume 5 | Article 153 | 11 



Korte et al. 



Desulfovibrio resistance to high nitrate 



using iTRAQ proteomics. Brief. Funct. Genomic. Proteomic. 5, 133-143. doi: 
10.1093/bfgp/eU025 

Riley, R., and Zachara, J. (1992). Chemical Contaminants on DOE Lands and 
Selection of Contaminant Mixtures for Subsurface Science Research DOE/ER- 
0547T. Washington, DC: US Department of Energy. 

Seitz, H. J., and Cypionka, H. (1986). Chemolithotrophic growth of Desulfovibrio 
desulfuricans with hydrogen coupled to ammonification of nitrate or nitrite. 
Arch. Microbiol 146, 63-67. doi: 10.1007/BF00690160 

Steinberg, C. E. W., Stiirzenbaum, S. R., and JVIenzel, R. (2008). Genes and envi- 
ronment - striking the fine balance between sophisticated biomonitoring and 
true functional environmental genomics. .SV(. Total Environ. 400, 142-161. doi: 
10.1016/j.scitotenv.2008.07.023 

Sunde, E., Thorstenson, T., Torsvik, T., Vaag, J. E., and Espedal, M. S. (1993). 
"Eield-related mathematical model to predict and reduce reservoir souring," 
in Proceedings of the 1993 SPE International Symposium on Oilfield Chemistry, 
(Society of Petroleum Engineers, Inc.), 449-456. 

Sunde, E., and Torsvik, T. (2005). "Microbial control of hydrogen sulfide produc- 
tion in oil resevoirs," in Petroleum Microbiology, eds B. Olliyier and M. Magot 
(Washington, DC: ASM Press), 201-213. 

Torres-Garcia, W, Zhang, W., Runger, G. C, Johnson, R. H., and Meldrum, 
D. R. (2009). Integrative analysis of transcriptomic and proteomic data of 
Desulfovibrio vulgaris: a non-linear model to predict abundance of undetected 
proteins. Bioin/ormfltics 25, 1905-1914. doi: 10.1093/bioinformatics/btp325 

Van Opijnen, T., Bodi, K. L., and CamiUi, A. (2009). Tn-seq: high-throughput par- 
allel sequencing for fitness and genetic interaction studies in microorganisms. 
Nat Methods 6, 767-772. doi: 10.1038/nmeth.l377 

Voordouw, G., Grigoryan, A. A., Lambo, A., Lin, S., Park, H. S., Jack, T. R., et al. 
(2009). Sulfide remediation by pulsed injection of nitrate into a low temper- 
ature Canadian heavy oil reservoir. Environ. Sci. Technol. 43, 9512—9518. doi: 
10.1021/es902211j 

Wahan, P. J., Allen, S., Shatsky M., Zeng, L., Szakal, E. D., Liu, H., et al. (2012). 
High- throughput Isolation and characterization of untagged membrane protein 



complexes: outer membrane complexes of Desulfovibrio vulgaris. J. Proteome 

Res. 11, 5720-5735. doi: 10.1021/pr300548d 
Wall, J., BUI Yen, H. C, and Drury, E. C. (2007). "Evaluation of stress response 

in sulphate-reducing bacteria through genome analysis," in Sulphate-Reducing 

Bacteria: Environmental and Engineered Systems, eds L. L. Barton and W. 

A. Hamilton (New York, NY: Cambridge University Press), 141-165. doi: 

10.1017/CBO9780511541490.005 
Wall, J. D., Rapp-Giles, B. J., and Rousset, M. (1993). Characterization of a 

small plasmid from Desulfovibrio desulfuricans and its use for shuttle vector 

construction./. Bacteriol. 175,4121-4128. 
Zane, G. M., Bill Yen, H. C, and Wall, J. D. (2010). Effect of the deletion of 

qmoABC and the promoter-distal gene encoding a hypothetical protein on 

sulfate reduction in Desulfovibrio vulgaris Hildenborough. Appl. Environ. 

Microbiol. 76, 5500-5509. doi: 10.1128/AEM.00691-10 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 27 January 2014; accepted: 21 March 2014; published online: 21 April 2014. 
Citation: Korte HL, Pels SR, Christensen OA, Price MN, Kuehl JV, Zane GM, 
Deutschbauer AM, Arkin AP and Wall JD (2014) Genetic basis for nitrate resistance 
in Desulfovibrio strains. Front. Microbiol. 5:153. doi: 10.3389/fmicb.2014.00153 
This article was submitted to Microbial Physiology and Metabolism, a section of the 
journal Frontiers in Microbiology. 

Copyright © 2014 Korte, Pels, Christensen, Price, Kuehl, Zane, Deutschbauer, Arkin 
and Wall. This is an open-access article distributed under the terms of the Creative 
Commons Attribution License (CC BY). The use, distribution or reproduction in other 
forums is permitted, provided the original author(s) or licensor are credited and that 
the original publication in this journal is cited, in accordance with accepted academic 
practice. No use, distribution or reproduction is permitted which does not comply with 
these terms. 



Frontiers in IVIicrobiology | Microbial Physiology and Metabolism 



April 2014 I Volume 5 | Article 153 | 12 



